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1. Introduction 


Research efforts on the development of the so-called large eddy interaction model 
(LEIM) have been devoted to the following: 

(a) setting up of dynamical equations for velocity and scalar quantities and 
correlations; 

(b) examining methods of solution of the dynamical equations with minimal 
modeling at the spectral level; and 

(c) examining data sets on mixing layers that are suitable for utilizing the LEIM 
as a metric for determining transport and production of turbulence quantities. 

The two-dimensional mixing layers under consideration fall under two categories: 

(i) the mixing fluids differ in velocity and density, the latter due to composition; and 

(ii) the mixing fluids differ in velocity and temperature. It is assumed that the fluids 
are incompressible and non-reactive. Thus, the dynamical equations for velocity and 
scalar fluctuations and correlations are common for density and temperature. 

The principal modeling required in the use of dynamical equations is for eddy- 

// 

eddy interactions. Alternative schemes of modeling are then of interest. Regarding 
the need for establishing a relation between velocity and temperature fluctuations, 
similarity may be invoked. In view of the nature of the LEIM, such similarity is most 
effective if invoked at the spectral level between velocity and scalar quantities. 

An attempt has been made at setting up a computational program for a 
developing mixing layer. There is considerable sensitivity to initial conditions in the 
prediction of turbulence quantities, although there is some success in predicting mean 
flow development. 


In respect of data sets, the two main problems are (i) the lack of adequate 
definition of initial conditions and (ii) the lack of data on mean flow development. 
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These preclude adequate comparison between predictions and experimental data. 

1.1. Background 

The large eddy interaction hypothesis (LEIH) is based on the work of J. Lumley 
(Refs. 1 and 2) on the rational representation of turbulence and its development. 
The LEIH rests on two postulates: (i) selected eddies are adequate for representation 
of a number of features of turbulence; and (ii) turbulence development in a flow is 
the result of the action of mean strain in the presence of eddy-eddy interactions that 
constitute turbulent transport processes and of viscous transport. 

The application of a model based on LEIH to curved wall boundary layers is 
discussed in Refs. 3-6. Other application to channel flows is discussed in Refs. 7-9. 
Recent developments on J. Lumley 's work are contained in Refs. 10-14. 

A general critique on experimental work on mixing layers can be found in Ref. 15. 

1.2. Outline of Report 

Chapter 2 of the Report deals with the rational representation of turbulence and 
its application to mixing layers. The method of obtaining predictions based on the 
LEIH model is presented in Chapter 3. An approach to the problem of predicting 
mixing layer development utilizing the LEIH is outlined in Chapter 4. 

2. THE LEIH AND MODEL DEVELOPMENT 

The LEIH is an outcome of J. Lumley’s proposal for rational representation of 
turbulence (Refs. 1 and 2), which has also been adapted for use by A. A. Townsend 
(Ref. 16) with his own formalism. A.A. Townsend had earlier (Ref. 17) proposed that 
turbulent shear flows have a "double structure" with a noticeable, well organized 
structure called large eddies and a less organized, smaller scale turbulence. The large 
eddies are said to contain only a small part of the total kinetic energy, but determine 
the turbulent-non-turbulent interface and hence the intcrinittency. There is expected 
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to be a large difference between the scale of large eddies and the scales of other 

eddies. The effect of the smaller eddies may be visualized as the presence of eddy 

( 

viscosity. A.A. Townsend was thus able to establish a scheme for determining the 
structure of large eddies. 

An objective definition for such large eddies has been provided by J. Lumley. The 
importance of such a definition is two-fold: (i) it leads to a model for inhomogeneous 
flows and (ii) it also provides a rational basis for the recovery of the characteristic, 
generally orderly structure of large eddies in various flows based on detailed space- 
time based measurements. In other terminology, it represents a means of examining 
time-wise unsteady, space-wise steady development of the structure of a given 
turbulent flow. 

2.1. Notation 

In order to make the notation compact and easily recognized, the following is 
utilized. 

Throughout Cartesian summation notation is employed, where repeated indices 
imply that the terms containing them must be summed over all possible coordinate 
indices. An overdot indicates partial derivative with respect to time and a dot inside 
the parenthesis (•) indicates that the quantity is a function of whatever are the 
relevant coordinates. 

The isotropic tensors are defined by 


5 jj = 1 and 0, if / = j and otherwise 
and 


tijk — 1 > — 1 and 0, when ijk form the sequence 123123, 321321 and otherwise. 
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Mean velocity, pressure and temperature are denoted by Uj, P and T, while 
turbulence velocity, pressure and temperature are denoted by uj, p and 6. 

All other symbols and operations are explained locally. 

2.2. Application of Proper Orthogonal Decomposition Theorem 

The turbulence problem is posed as one of identifying a structure in a random 
vector field. The proper orthogonal decomposition theorem of M. Loeve (Ref. 18) is 
invoked for application. 

The vector field Uj is expanded in terms of a candidate structural quantity, 
<f>f , n = 1,2,. ..etc., and coefficient, a n , as follows : 

(2.1) 
( 2 . 2 ) 


«,<•> = E «„ *,"(•> 

where, denoting the complex conjugate by ( )*, 
ar„ = f«i(*) <£"’(•) <*(•) 


and the random coefficients are uncorrelated: 




0 , n * m 
\ (n) , n = m 


(2.3) 


It is claimed as a hypothesis that is occurring in a recognizable form in a given 
ensemble of random vector fields, Uj. One considers the quantity, namely 

U;(.) «,(.) d(.) 

« = T7 (2.4: 

<M*) */<•> d(OJ 


and its statistics. The best <j\ is then taken to be that which gives the largest 
magnitude of a in some average sense, the sign being irrelevant. In other words, </> j is 
selected such that | a | 2 is maximized. According to the calculus of variations, the 
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expansion of (2.1) is then optimal since, on account of | ex | 2 being maximized, each 
term leaves as little as possible of Uj(') in the succeeding terms. Noting that the large 
eddies are physically well "separated" from the rest, the first term of the expansion 
can be claimed to be the representation of large eddies. 

Defining Rjj as 

(2.5) 

one can write 

*//(•. •') =EX'"V/" > (.) */">(•') (2.6) 

It may be noted that in (2.3) are non-negative and their sum is finite. 

Furthermore, 0( n ) are ortho-normal, i.e., 

(^(•J d(>) =b pq (2.7) 

This is to say that the structures of various orders have nothing in common with one 
another. 

Now, one considers the possibility of Uj or Rjj being homogeneous in certain 
coordinate directions and inhomogeneous in the remaining. For example, the vector 
field may be homogeneous in x 3 and inhomogeneous in Xi and x 2 , as in a two- 
dimensional boundary layer or a two dimensional mixing layer. Then one can invoke 
the Harmonic Orthogonal Decomposition Theorem (M. Loeve, Ref. 18) and note that 
expansion in harmonic functions is possible in any homogeneous directions, while the 
proper orthogonal decomposition applies in the remaining. 


It is clear that there is a difficulty in a totally homogeneous flowfield. Considering 
a wave number-frequency space, the more broadband in wave number and frequency 
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is the structure, the less representative will the definition be. On the other hand, if 
the flowfield is totally inhomogeneous, one has to appeal to </>[*) . In the case that 
there is at least one direction in which the field is homogeneous, the flow can be 
represented in terms of a denumerable set of motions. The large eddy then is that 
organized motion of largest scale in the inhomogeneous directions and possessing the 
greatest energy. A strong peak is expected at a wave number-frequency combination 
away from the origin and that corresponds to the large eddy. 

2.3. Some Comments 

2.3.1. The eddies in the turbulent flowfield are represented in terms of flow quantities 
treated as random vector fields and hydrodynamic/thermodynamic quantities, such 
as pressure, temperature and concentration, as random fields. Various correlations 
between velocities and between velocities and scalars can then be constructed. 

Each field is represented in terms of a structure and associated random 
coefficients representing amplitudes. Upon decomposition, a flowfield with some 
inhomogeneities can be represented in terms of a denumerable set of eddies, each of 
which is rational being optimal in the expansion. Each contains the maximum energy 
it can and each is of the largest scale possible. There is no similarity among them. 
The largest of them is the large eddy. 

2.3.2. As stated earlier, the Orthogonal Decomposition Theorem reduces to the 
Harmonic Decomposition Theorem when the flowfield is homogeneous, and also, 
when it is stationary or periodic. This implies that in any direction in which the 
flowfield is homogeneous, stationary or periodic, the flowfield in that direction can be 
expressed by Fourier analysis. In the remaining directions, one has to use the Proper 
Orthogonal Decomposition Theorem. 

J. Lumley points out (Ref. 2) that eigen functions are not confined to one region 
of space and time and therefore they cannot correspond to physical eddies as they are 
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usually visualized. 

However, he also points out (Refs. 1 and 2) that, following a suggestion of A. 
Townsend, one can utilize an extension of the shot-noise-efTect expansion of S.O. Rice 
(Ref. 19) and obtain an unambiguous representation of eddies. Now, any 
homogeneous function may be written in the form, namely 

«,(■*) = f fi(x-x) g( x) dx’ (2.8) 

where f; is a deterministic function and g is a stochastic function that is white, i.e. 
uncorrelated in non-overlapping intervals. Thus, the flow is decomposed into 
characteristic eddies, gj, occurring at uncorrelated intervals. One can then use the 
extended Campbell Theorem and write 


Rij(Z) =\fMfj(x+$) dx (2.9) 

Taking the Fourier transform, one can express the spectrum as 

* ij (k)=fif j (2.10) 

where ( A ) indicates the transform. 

2.3.3. L. Sirovich (Ref. 20) points out that another method, namely the method of 
snapshots or strobes may be utilized in place of the shot-noise method. We shall 
return to the alternative method later. 

3. DYNAMICS OF LARGE EDDIES 

In any turbulent shear flow it is hypothesized that the large eddies interact with 
themselves and with all the other eddies in the presence of mean strain in the fluid 
flow. A dynamical equation can be constructed for any parameter associated with the 
large eddies, wherein a balance is invoked among production, advection, diffusion 
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and dissipation processes. Turbulence development in the flow, including mean flow 
development, is assumed to be governed by the dynamical equation along with other 
relevant flow describing equations. 

As noted by J. Lumley (Ref. 1), for the completely inhomogeneous case, 

«>,(•) = X (m) 4>, (m) (*) (3.1) 

which operation is a consequence of Eqs. 2.1 - 2.3. Substituting Eq. 3.1 into the 
Navier-Stokes equations, one obtains a dynamical equation as follows. 


4>! n) + U u 4>/ n) + W Uj + J] (4? ap «,) /a 1 ‘ n) + £ K* /\ (B > • */«> = - - tt/"* + .*$ (3.2) 






^7 } = o 


(3.3) 


= * (m) (3.4) 

Here Uj, p and 7r represent the mean velocity, density and pressure, and n and m 
represent the order. The first two terms on the left side of Eq. 3.2 represent the 
advective derivative. The third term on the left is the result of interaction with the 
mean strain. The next term denotes interaction with all other eddies. The two terms 
on the right hand side of Eq. 3.2 are related to the pressure gradient and the viscous 
action. 
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Now, the central problem of Eq. 3.2 is again the so-called closure problem because 
of the term denoting eddy-eddy interactions, which requires some form of 
approximation. Several alternatives are as follows, although others may be feasible: 

(i) a direct interaction approximation of the type given by R.H. Kraichnan (Ref. 

21); 

(ii) a random phase approximation of the type given by S.F. Edwards (Ref. 22); 

(iii) a type of approximation given by W. Heisenberg (Ref. 23) which, although 
proven unsatisfactory in the case of homogeneous turbulence, has shown some 
promise in inhomogeneous flows; thus, A.S. Monin and A.M. Yaglom (Ref. 24) 
indicate that the higher order eigenfunctions may be considered as acting like a 
viscosity on the lower order ones; the higher and the lower order terms are 
assumed to be sufficiently different from each other based on their spectral 
separation. This has been pursued extensively in recent times by N. Aubry 
(Refs. 10-12); and 

(iv) direct modeling of the term on the basis of its relation to turbulence transport 
(Refs. 1,2,3 and 7). 

It may be recalled again that L. Sirovich follows a slightly different procedure. 

3.1. Partially Inhomogeneous Case 

In shear flows where there exist approximately homogeneous directions, for 
example Xj and x 3 directions while x 2 is inhomogeneous, and where the turbulence is 
stationary, one can write for the fluctuating field, as stated earlier, the following; 
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u, = ^ exp | 'nkjXj + wl) • 4>, in) a„ dk\ dk 3 du 
n 

= j fi gi d* 


(3.5) 

(2.8) 


and 

f=\ti l) dk t ... ( 3 - 6 ) 

It will be observed that Eq. 3.6 denotes the deterministic function. The dynamical 
equation for i// 1 ) then becomes the following: 

+ tf'/J 1 Uj + Ujj \pf ]) - - p ( P + (^*) ) UipiDt (3.7) 

where T s is a time scale, proportional to the ratio of turbulent energy to dissipation 
rate. It should be pointed out that viscous dissipation is not included in Eqn. 3.7. 
Equation 3.7 has the same features as Eq. 3.2. 

3.2. Eddy-Eddy Interactions 

It is clear from Eq. 3.2 that in incompressible flow, whether or not scalar 
quantities are present such as temperature and species concentration, the equations 
are closed insofar as pressure and viscous forces are concerned. Of those pressure can 
be established by including Poisson’s equation. The only term that is unclosed is that 
pertaining to eddy-eddy interactions. 

First, it is useful to note that Eq. 3.2 is a coupled, infinite set of eigen value 
equations. In principle, they should yield the eigenvalues a*a p a q /\( D ) when solved 
simultaneously. The solutions yield moments of higher order in terms of those of 


lower order. 
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A brief discussion follows on the modeling of the interaction term. 

3.2.1. Lumley’s original suggestion 

Based on W. Heisenberg’s suggestion (Ref. 23), J. Lumley introduced the 
hypothesis that a p </>[ p ) extracts, through some form of viscous action, energy from 
n<p, as a function of time at appropriate scales. One can then write: 

@q 0 , p 5 ^ q 

and 


— o' 


Op Oq 




= 


®np &ij 


+ \ {n) v { T p) 


fj. (n) 
( ViJ 




p,q > n 


(3.8) 


Unfortunately, there is no mechanism for the extraction of energy in the formalism 
itself. Therefore, at best one can consider the first mode and write the following. 


+ U,j <t>; 


<i) 


+ Uj 


a*a\ 


j ( 1 ) 


A d) 

<t>j = 



+ (-+ E *H’ 


(3.9) 


and 




One introduces in addition the Poisson equation. 

Three important remarks of J. Lumley regarding the foregoing are as follows. 
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(i) The interaction of close wave numbers does not seem to be fully accounted 
for. It is of course possible to proceed assuming that there is no correlation between 
even close wave numbers. 

(ii) The amplitude parameter, namely (a^aiJ/Xj cannot be determined. The 
argument is that amplitude may change without a change in structure while the 
disturbance remains neutral at the same However, since amplitude change does 
cause structural change, one has to adopt some other argument, such as the necessity 
of a maximum value of t*r, to determine the amplitude factor. 

(iii) One could write 


= 0 ^ ( 310 ) 

where (3 is in the nature of a structural constant, a skewness factor. 

Equation 3.10, along with a value for i/j, provides a basis for examining shear 
flows. 


3.2.2. An eddy viscosity model 

The nonlinear eddy interaction term, representing turbulent transport and 
dissipation, may in general be expressed as a function of the local value for the 
large eddy and of its first few derivatives (Ref. 25): 


£ ♦,<*> =>(*,<", */". <C) 

i 


(3.11) 


Various types of approximations to transport can be obtained as special cases of Eq. 
3.11. 

In particular the commonly-employed gradient diffusion hypothesis (Ref. 26) 


yields the following. 
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oo 

L i(p) j(q) , (I) 

<t>j = *ik <t>kj ( 3.12 

P<i= 1 

where ejj is the anistropic eddy viscosity. 

Since the left hand side of Eq. 3.12 is symmetric with respect to i and j, the eddy 
viscosity can in the most general case be a fourth order tensor. However even the 
fourth order viscosity, e^i, is known to be not satisfactory in all respects, such as the 
violation of the condition of invariance with respect to rotation of the coordinate 
system (Ref. 27). Usually a second order eddy viscosity, fjj, is assumed to be sufficient 
to account for the anisotropic nature of an inhomogeneous turbulent flow. 

It may be pointed out that the use of eddy viscosity implies an infinite velocity of 
transport for the transfer of energy from large eddy interactions to subsequent 
smaller eddy interactions. ' 

The application of this approximation to curved wall boundary layers is discussed 
in Refs. 3 and 4. 

3.2.3. Finite velocity of transport 

An alternative to the (generally employed but questioned) gradient diffusion 
model is to assume that some form of "bulk convection" forms the basis for turbulent 
transport (Refs. 17, 28 and 29). 

Introducing anisotropic, transport velocity scales, Vjj, to represent transfer of 
energy from a high turbulence intensity region to one of lower turbulence intensity, 
one can write as follows: 


£ « V, 

P’l 7 = 1 


( 3 . 13 ) 


The transport velocity scales must again be symmetric. 
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The application of this approximation to curved wall boundary layers is discussed 
in Refs. 3 and 5. 

It may be pointed out that the foregoing method of relating the nonlinear term to 
local turbulence intensity through the use of a finite transport scale converts the 
dynamical equation into a hyperbolic one. Refs. 3 and 5 then discuss a characteristic 
field velocity associated with the characteristics. Considering other assumptions 
introduced in Ref. 29. one can write the characteristic velocity in the following form: 


V* = i (v 2 c -2a, (3.14) 

where V c and aj are empirical quantities introduced when modeling diffusion of 
turbulence kinetic energy and turbulence energy itself, respectively. The velocity V* 
then represents the characteristic speed associated with local changes in shear stress. 

It is interesting in this connection that one can also deduce a characteristic 
velocity for turbulent transport of the temperature field in a mixing layer based on 
Ref. 30. The characteristic velocity becomes 


V* = v 2 


(3.15) 


where v is the constant of proportionality introduced in modeling triple correlations. 
It is also equal to the turbulent normal stress. Then, assuming isotropy in transport 
velocity scales, one can write 


V* 



(3.16) 


The role of the transport velocity scales, Vjj, in the transport of eddy-eddy 

_2 

interactions is analogous to that of v in the transport of temperature fluctuations, 
but the two velocity scales are different. 
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One can also examine the transport of velocity-temperature correlations. 

The precise analogy between characteristic velocities for momentum and scalar 
quantities such as temperature is yet to be explored. We will return to this subject 
later. 


3.2.4. Aubry’s adaptation 

The main thesis of N. Aubry (Ref. 10) is that energy transfer from mode to mode 
can be modeled by a representation adapted from Heisenberg’s spectral model (Ref. 
23). Eventually, this model is utilized for establishing the dynamical processes 
associated with the characteristic structures of turbulent shear flows; we return to 
this topic in Section 3.3. 

Heisenberg's model in isotropic and homogeneous turbulence is based on the 
following assumption: in a given energy spectrum of turbulence, the energy transfer 
to wave numbers larger than a cut-off value k may be set equal to an effective 
transport coefficient ^(k, T), where T(k) is a nonlinear spectral transfer term; the 
transport coefficient is a measure of the action on eddies of wave number smaller 
than k. There is a notional similarity with augmentation of thermal agitation. 

Aubry then introduces the assumption that one can average the influence of the 
small scales on the large ones; to quote, "The smaller scales extract energy from the 
larger ones by a global viscous action on each unresolved mode." 

To proceed, one may adopt the same notation as N. Aubry: contributions less and 
greater than the cut-off between the large and the small scales are denoted by 
subscripts < and > respectively; an average value is indicated by < >; and 
additional subscripts to the average denote ranges over which averages are obtained. 

Assuming that the small scale stress tensor is again linearly proportional to the 
strain rate tensor of the large scales, one can write the following. 
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T ij> = ~2v t S iJ< 


(3.17) 


where 

Sij < = ~ j («/<,; + «/<,/) ( 3 - lg ) 

In terms of mean and fluctuating components of velocity, U and Uj, for instance, 
one can set up a dynamical equation as follows, exactly as was done in Eq. 3.2: 


«,< + u iKj Uj + U i< j*u j< + u i<J uj - <u i<J u j< > 

= ~ “ P + ^ &ij (<“k> u k> > ~ «“k> »k> ») + (" + i't) u,<jj (3- 19) 
p 3 x J 

which applies to the small wave number part. 

There arise two terms in Eq. 3.19 that require some form of modeling. One of 
them is combined with the pressure term. In order to calculate this, it is assumed 
that the fluctuations of the kinetic energy of small scales are proportional to the rate 
of loss of energy by the large scales to the small scales. 

The other term requiring modeling in Eq. 3.19 is associated with v <p. Again, it is 
assumed that can be made proportional to the kinetic energy distribution of the 
higher modes over the shear flow under consideration and a characteristic length 
scale; the latter is set proportional to the ratio of kinetic energy and dissipation 
energy. 

It is clear that in both cases, the assumptions are the usual ones and therefore 
appear as artifices that are justified by their success. A further controversy pointed 
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out by J. Lumely relates to the neglect of any influence of the mean strain in Eq. 
3.17. Finally, it is useful to note that in sub-grid scale modeling, P. Moin and J. Kim 
(Ref. 31) introduce the following expression for the small scale stress tensor. 


T U> = “ v t (Sy ~ <Sjj >< ) - v' T < Sy>< (3.20) 

which can be compared with Eq. 3.17. 

3.2.5. Developments in the use of (S, i/p). 

It will be recalled that in Section 3.2.1, a suggestion was made that the amplitude 
parameter can be related to a factor, /?, Eq. 3.10, which is in the nature of a skewness 
factor. Neither (3 nor Uj are strictly "physical" quantities and we introduce the 
notation S and z>r to represent them as in the nature of the associated quantities. 

S.K. Hong (Refs. 7, 8 and 9) has examined the possibility of simulating channel 
flows with a model based on the LEIH utilizing S and t/j as the transport and 
dissipation factors. 

The formalism may be summarized as follows. 

The mean velocity is decomposed into an ensemble average component and a 

fluctuating component. The latter is written in terms of orthogonal functions, 4>/"\ 
n = 1, 2, 3 Thus 

OO 

«,(*.') = Yt (3.21) 

n = 1 

<t>i q) dx dt = 0, p * q, i = 1,2, 3. (3.22) 

where c* n are random coefficients with units of velocity and are uncorrelated with one 
another. Ilence, one can write the following. 
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s: = 0 


and 


«A = ^ (n) K, 

where X( n ) are all positive. Now, the </>(“) are assumed to be orthonormal functions, 
implying 

\<t>, iP) <t>l q) dx dt = 8 pq 

and, further, that none of the <f >\ D ) is identically zero. It can then be shown that a„ 
are Fourier coefficients given by 

a„ = j Ui(x.t) dx dt (3.23) 


where i *= 1. 

If one considers the integral of the turbulent kinetic energy of the Uj —component 
across the whole flow, one can write 


- u 2 </*<* = - y x (n) 

1 l LJ 


' n= 1 


(3.24) 


The n— th mode distributes the energy in space and time according to its functional 
form. 

In the case of turbulent flow, a dynamical equation corresponding to the Navier- 
Stokes equation becomes the following: 


«; + Uj Ujj + Ujj uj + (u.Uj - UjUj),j = -~p,i + p Ujjj 


(3.25) 
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Introducing the orthogonal decomposition of Eq. 3.22, after some algebraic 


manipulation, 


one obtains a dynamical equation for 



as follows; this is similar to 


Eq. 3.8: 




+ Uj $ 


,0 + U'J $$ 


+ 




(3.26) 


where 


$i n) = vx^ <t>: 


(»> 


(3.27) 


( L-f£ E 


p=I <?=1 


QQ. <*P a q 



(3.28) 


and 

f <n) = -p- 1 (3.29) 


It may be observed that the pressure fluctuation term, 7r^, has no restrictions other 
than detenninacy; in fact, one can utilize the Poisson equation simultaneously with 
Eq. 3.26. 


For the first mode, the dynamical equations for and 7T^ become the 


following. 


r-d) 

<p, 


+ Uj Uij <£ (, V 



+ V 0 


' 0 ) 

iJJ 


(3.30) 
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and 




Again, 


( 3 . 31 ) 


}-E 

p - 1 


y a ' a p^ 
[x^X^X^]* 


*} p) f* w 


and 


( 3 . 32 ) 


a, 


00 00 

f 1 = V V 

1 J <,=, x ( 1 >x<'>x<«> 


f/«’ 0* (p> 


( 3 . 33 ) 


The central problem, then, is to model the term in Eq. 3.32 and Eq. 3.33. Either 
term represents modal interactions giving rise to transport and dissipation. It is 
proposed to write 


( )-■ 


«? 


x < 1)3/2 


r-o) .re) a (') 

<Pi <Pj ~ V T <Pij 


= S </> / 1 * _ V T fy) 


r<>) 


( 3 . 34 ) 

( 3 . 35 ) 


wherein S and i/p are in the nature of skewness and eddy viscosity. Considering S as 
a structural factor involved with transport, is? i9 associated with damping. In general 
S is close to one numerically except for the influence of inhomogeneities in flow. It 
will be observed that the skewness is associated with the nonlinear term in the 
dynamical equation. The term itself, therefore, must not be neglected. Once the term 
is retained, it becomes essential to have an additional dissipation in the form of //-p- 
The method, accordingly, requires two parameters to be selected. There is no direct, 
rational procedure for defining or selecting the parameters. However, an indirect 
justification can be obtained and is discussed later. 
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Now, considering a shear flow, such as a channel flow in which the inhomogeneity 
may be considered to be confined to the y-direction, in standard Cartesian 
coordinates, one can define spectral functions as follows. 


(k h y,k 3 ,t) = Jj[^ ^ <U I** 1 ,t) exp | -i (£|jc + fc 3 z)jj dx dz 


(3.36) 


and 


TT (k t ,y,k } ,t) = — — 5 || Vx 7T (I) exp | -f(k { x + £ 3 z)j] dx dz (3.37) 

( 2w ) 

A similar approach can be adopted for a two-dimensional mixing layer, which is 
discussed further in Section 5. 

4. RECOVERY OF CHARACTERISTIC STRUCTURES 

The existence of some form of structure in turbulence has been postulated for 
several decades based on various types of observations. However, studies of their roles 
in different flows and of their physical structure has become increasingly fashionable 
over the past fifteen years or so, as pointed out by H. Liepmann and expounded in 
Refs. 32. There are obviously controversies in regard to the importance of their role, 
their existence in fully-developed turbulence and in various types of shear flows and 
their strength. Nevertheless in actual fluid machinery, where turbulence may not 
attain a fully developed state, and in the case of mixing layers in particular, there is 
reason to believe that (Kelvin-Helmholtz) instability-generated large structures may 
display coherence and a specific role based on their structure and the gradual break- 
down of their structure. Regarding fluid flow machines, it is however, necessary to 
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point-out that various (inherent) disturbances may in fact disorganize the structure. 
The mixing layer has, in any case, been postulated by A. Roshko (Ref. 33) to be quite 
distinctive in this respect among inhomogeneous shear flows. 

J. Lumley and F. R. Payne (Ref. 34) obtained the large eddy structure in the 
wake of a cylinder within the premise of the theory put forward by J. Lumley, Ref. 

I. It is the first accomplishment of the demonstration of (the large) eddy structure 
starting from classical, statistical, structure theory of turbulence and utilizing the 
most rationally-based representation of the structure. It is important to point out 
that no phase information was generated in that early work and therefore nothing 
was said about the coherence or the phase-relations governing the structures in the 
mean fluid flow. In any case, the possibilities for determining a structure became 
clear and A. Townsend has a critique, in proper perspective, of this approach in Ref. 
16. 

A method of extending the earlier analysis to the determination of coherent 
structures has been outlined in Ref. 3. A modification to the shot-noise approach of 

J. Lumely utilizing a snap-shot approach has been outlined by L. Sirovich in Ref. 20. 
A modification to the recommendation in Ref. 3 has been worked out by N. Aubry, 
Refs. 10 and 11. One should also note the important investigations reported by P. 
Moin in Ref. 35. Finally, M. N. Glauser and W. K. George (Ref. 36) have attempted 
a recovery of structures in an axisymmetric jet, utilizing experimental results along 
with orthogonal decomposition analysis. It is considered important to review this 
body of investigations if only briefly. 

4.1. Lumley’s suggestion for recovery of phase information 

In Section 2.3.2, the spectrum of <£jj has been given by Eq. 2.10. That equation 
permits the determination of f to within a phase angle given by the following. 


23 


fa = <t>aa ex P *®(*) 


(4.1) 


In other words the phase or the overlap and spacing information cannot be obtained 
from second order statistics. It is then shown how the phase information can be 
obtained from third order statistics, 

u(t) u{t + 7 ,) u(t + T 2 ) 

The procedure is elaborated in Ref. 3. 

In order to obtain overlap and spacing, one again has to utilize higher order 
correlations of g,(t). Once the phase of f has been determined, one can Fourier 
transform u,(t)and the individual values of the transform can lead to obtaining the 
transform of g(t). 

One point to note is that an important piece of information is required, from 
experiments or empiricism, about g, the stochastic function. 

4.2. P. Moin’s Application to Channel Flow 

The objective (Ref. 35) was to extract coherent structure details for a channel 
flow in a two-dimensional plane (x, y), y being the inhomogeneous direction normal 
to the wall. The structures were assumed to be sprinkled randomly in the spanwise 
direction. No attempt was made to obtain phase information. The numerical results 
were compared with those from large eddy simulation model. 

It is not clear how the mean flow is prescribed or obtained in the foregoing 


calculation. 
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4.3. Sirovich Method 

In Section 2.3.3. we mentioned briefly the recommendation of L. Sirovich for the 
use of the so-called snapshot method (Ref. 20) as an alternative to the use of a shot- 
noise method utilized by J. Lumely. 

The methodology depends upon the availability of sufficiently extensive 
information from numerical or physical experiments in which highly resolved flows 
are followed in time. A coherent structure is associated with the eigenfunctions of the 
correlation operator. The main problem then is establishing a method for deriving 
the eigenfunctions. Several aspects of the method become clear in the next section. 

It may be stated that no numerical results are presented by L. Sirovich. 

4.4. N. Aubry’s Approach 

We continue here the discussion presented in Section 3.2.4. 

D. Ruelle and F. Takens (Ref. 37) proposed a model of turbulence as a 
deterministic, chaotic regime of flow reached after a small number of bifurcations, 
rather than one involving indefinitely many modes. The proposal was based on a 
possible connection between dynamical systems theory and turbulence. It is generally 
accepted that chaotic, dissipative, dynamical systems eventually display features of a 
strange attractor that is of low dimension. L. Sirovich suggests an optimal 
description for the attractor and also some form of an upper bound to the dimension 
of the attractor. Now, N. Aubry (working quite independently) has investigated the 

question of "whether the dynamics of a complex turbulent flow can be described 

by a finite, possibly small, number of modes". The complex turbulent flow chosen for 
illustration is the wall region of a turbulent boundary layer in a channel flow, where 
large scale, horseshoe shaped, structures have been shown to play an important role. 
These organized, coherent structures appear in a fine-grained turbulent background. 
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A.K.M.F. Hussain (Ref. 38) has stated that "the motion of coherent structures .... is 
likely to be low-dimensional". The problem is to extract deterministic structures in a 
description of turbulence based on the rational representation of J. Lumley. 

One question, raised by L. Sirovich, pertains to the fact that coherent structures 
evolve in time. The original suggestion of J. Lumely (Ref. 2) was to construct a 
coherent structure as an appropriate superposition of eigenfunctions, albeit utilizing 
higher-order statistics for calculating the coefficients. This approach is unsuitable for 
determining evolving structures. L. Sirovich therefore utilizes Galerkin projection, in 
addition, as stated earlier, to abandoning the shot-noise method. N. Aubry also has 
utilized Galerkin projection along with proper orthogonal composition. The Galerkin 
projection is an effective means of reducing truncation errors and also, leads to a set 
of ordinary differential equations for the coefficients. 

A significant point here is the claim to establishing "a unique relationship between 
the correlation tensor and the unsteady flow that produces it," which had been 
thought impossible by B. J. Cantwell (Ref. 39). 

The method is implemented by N. Aubry to determine the characteristic eddies in 
three dimensions together with their temporal dynamics. The flow is considered 
steady, but periodic in two directions, while being inhomogeneous normal to the wall. 

An issue in the method is the manner of prescribing the mean flow. A mean 
velocity profile is obtained from a physical experiment for a fully developed channel 
flow. However, this cannot be introduced into the formalism wherein one obtains, 
starting from the Navier-Stokes equations, the following for a fully-developed 
channel flow: 


<u x u 2 > 2 = - P ] p, i + vU, 22 


( 4 . 2 ) 
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and 

<« 2>.2 = “ P ~ l P ,2 ( 4 ‘ 3 ) 

Expressing the wall mean pressure gradient in terms of friction velocity and 
integrating Eq. 4.2., one obtains 

< 44 > 

where H is the half-height of the channel and u r , the friction velocity. N. Aubry 
points out that the U-profile introduced needs to be consistent with Eq. 4.3. In other 
words the velocity profile must show the interaction between the Reynolds stress and 
itself. Physically, the shear stress term is related to the intensity of the structures, 
which is inversely proportional to Reynolds stress. However, this raises a problem in 
selecting the U-profile objectively. 

It may be pointed out that both P. Moin and N. Aubry obtain the structure 
based on the large eddy represented by the first model. The question of choosing and 
using a certain number of modes is addressed in the next section. 

4.5. Multiple- mode Analysis 

J. Lumely’s rational representation of turbulence can be summed up as follows: 
the large eddy is that structure which has the largest mean square projection on the 
velocity field. It can be defined through orthogonal decomposition of the field of 
interest. Its structure can be established utilizing a dynamical equation. In defining 
it, one can introduce a truncation at any low mode number; for example, the first 
mode can be considered as the large eddy. On the other hand, one may want to 
include several modes. Each mode must contain the largest amount of energy 
possible. However, there is no rational basis for determining the energy transfer from 
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mode to higher mode. In principle, one can consider a set of dynamical equations, a 

set of eigenvalue equations, which is solved to obtain eigenvalues of a n a p a q /X^. In 
practice, it is simpler to obtain the required information from physical experiment. 
This procedure of utilizing data from physical experiments has been utilized by W. 
K. George (Refs. 40-42) for determining coherent structures in an axisymmetric jet 
mixing layer. No attempt was made to obtain phase information, and therefore, only 
the local structure has been obtained. 

In the axisymmetric jet studied there is stationarity with respect to time and also 
periodicity and symmetry in the azimuthal direction; inhomogeneity exists in one 
direction. In the inhomogeneous direction one uses orthogonal decomposition after 
fitting the homogeneous directions with harmonic eigenfunctions. 

A scalar two-dimensional version of orthogonal decomposition has been derived: 


\<t>(x,x',o)) <j> n (x',ui) dx' = X (ra) (to) 0 <n) ( x,co) 


(4.5) 


The characteristic eddy is determined from the dominant eigenfunctions as 





0 (I) (JT,CO) 


rfco 


(4.6) 


Experimental data are utilized for the cross spectrum. 

Several modes have been superposed to examine the resulting structure. That 
combination which appears most acceptable compared to experimental data is 
retained. It is clear that no rational procedure is available for determining the higher 
order modes or selecting a combination of them. 
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We have referred to this fact earlier in Section 3.2.4. Referring to Eq. 3.8, 
assuming from experimental observations that <f ) is weak, a solution can be 
constructed in the form: 


lO) aS^ -L * c 

>, <Po,i + t<P hi + 6 


2 , 0 ) , 

+ 


(4.6) 


where 

e =ar«?/x <l) 

In Eq. 4.6, the first term determines the (neutral) stability of the mean motion to 
small disturbances. The parameter i/j can be selected on this basis (Ref. 16). 
However it turns out that at least the second term must be retailed when the flow is 
such that the mean velocity profile is stable with respect to small disturbances. 

Now, a similar optimization can be carried out, in principle, to determine a 
succession of modes; however, as stated earlier, this may not be effective. 

The evolution of modes within the framework of Navier-Stokes equations can be 
determined noting that Fourier decomposition of the equations leads to a triad 
structure in Fourier space (Ref. 43). The triad interaction represents the nonlinear 
process of energy exchange between the modes. The interaction depends upon the 
triad geometry, the orientations of the amplitude vectors with respect to the triad 
plane and the phase relations among the modes. A basic study on the generation of 
various modes and determination of their energy content is reported in Ref. 44. 

It may be added here that W.K. George and N. Aubry also construct the 
streamline pattern within a coherent structure of the flows considered by them. 

Finally, one may note that several scales of interest can be constructed utilizing 
the values of . 
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The scale obtainable most directly is the average Taylor scale (Ref. 1), 



Based on spatial velocity autocorrelation, one can define the Taylor microscale in the 
form 



( 4 . 8 ) 


The integral length scale can be written as 

L e = fiY 2 je ( 4 - 9 ) 

where e is the dissipation rate of energy. A time scale can then be defined in the form 


Lg/Uj 7 l0 


5. A MODEL FOR A MIXING LAYER 

We consider a mixing layer, as shown in Fig. 1, involving two-dimensional mixing 
of two streams of Newtonian, compressible, turbulent fluids differing in mean velocity 
and temperature. The general interest is in predicting (a) rate of spreading of the 
mixing layer, (b) intermittency at the two boundaries, (c) coherent structures with 
their streamline pattern and (d) critical interaction between coherent structures. Our 
objective here is limited to examining a formalism of LEIH that is applicable to a 
mixing layer. 

5.1. Describing Equations 


Two sets of describing equations for the mean flow quantities are provided in 
Appendix I: (1) mass-averaged compressible flow equations (Ref. 45), and (2) 
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equations with no change in mean density and neglect of density fluctuations as small 
(Ref. 16). The latter set represents a physical case of mixing of two flows with a small 
difference in temperature or concentration, and therefore of interest in addition to 
being a useful test case for the application of the proposed model. 

It may be pointed out that it has repeatedly been stated, for example in Ref. 47, 
that density non-uniformity influences only the mean flow quantities and turbulence 
is unaffected by density fluctuations in thin mixing layers. These assumptions are 
equivalent to a local equilibrium approximation and can be applied, with great 
computational convenience, whereby they are shown to be adequate. 

5.2. Representation of turbulence 

We consider two main types of representations: (1) as in earlier work, large eddies 
in one group and all other eddies in another group, and (2) large eddies, energy 
containing eddies and all other eddies. The latter has some resemblance to the model 
of Ref. 46. 

It is also possible to visualize yet another representation with active, inactive and 
all other eddies. This is not of particular interest in the case of mixing layers. 

The two representations of interest in mixing layers are discussed in Appendix II. 

5.3. Temperature Field 

In the given formulation, it becomes necessary to model the eddy-eddy 
interactions separately with respect to velocity and temperature. Including turbulent 
viscosity and heat diffusivity, this entails assigning three skewness-like factors for 
velocity, temperature and correlations of the two. The Poisson equation for pressure 
fluctuation is unaltered in incompressible flow. 

Regarding the skewness factors, one possibility is to invoke a type of structural 
similarity between velocity and temperature that has been discussed extensively in 
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Refs. 48-49. In summary, the analogy suggests that 


q 2 /e 2 


dT 

dy 


dU 

dy 


is a constant across different parts of the mixing layers provided the averaging is 
done over appropriate time scales. 

The possibility of using such an analogy thus rests on investigating eddy-eddy 
interactions over a number of time scales starting with the time scale associated with 
the generation of large scale structures, through the time scale based on a 
characteristic length scale and a velocity of motion of large scale structures, 
ultimately to dissipation time scale. An approach of this type is discussed in Ref. 50. 
A method of incorporating such an approach into the large eddy interaction model 
was under investigation at the close of the grant period. 
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APPENDIX I 
BASIC EQUATIONS 

1. General Describing Equations 

Standard terminology as in Reference 45 is employed throughout. 

Considering a multi-component gas mixture involving N species and M elements, 
the equation for conservation of individual species may be written as follows. 


d d d 

- (py>) + — (p = — 

Ot OX k OXfc 


pD 


dxi 


+ w} i — l, 2, .... N 


0 . 1 ) 


It will be observed that a single diffusion coefficient is utilized. The last term in the 
foregoing equation represents the chemical source term. 

The equation for overall conservation of mass becom s: 

^ + -J- ( pUk ) = 0 
at dx k 


The conservation of momentum is given by the following: 


d d dp d 

T ( P «<) + T- (p «*«,)=- — + — T ik 
dt dXfc oXj ox^ 


(1.3) 


here p is the hydrostatic pressure, a thermodynamic variable. The viscous tensor is 
given by 


dxj dxjj 




Tij = P 


(1.4) 
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The conservation of energy can be expressed as follows: 


d d dp d 

- ( ph ) + — (pu k h) = — + — 
dt dx k dt dx k 



The mixture enthalpy, h, can be expressed by the relation 


h = £ y, hi 

i= 1 

where 


( 15 ) 


( 16 ) 


hi =c p T +4 ( 1 . 7 ) 

which accounts for the variations over the range of interest in temperature and heat 
content changes due to chemical transformation. The Prandtl and the Schmidt 
numbers are denoted by c and Sc, respectively. 

In most flow of gases involving chemical reaction, one can expect a density 
change, the density variations arising due to inhomogeneities in temperature or 
composition. It is then common practice to utilize Favre averaging, wherein all 
quantities, except pressure, are mass-averaged. The equations based on Favre- 
averaging corresponding to equations 1.1 - 1.3 and 1.5 become the following in the 
absence of variation with time. 


a a 

— ( pu k Y i ) = — 
OXf. dx k 


( dY, 

(pD- pu k 1 

\ dx k 


Y 7 + 


(1.8) 


~ (p u k ) = 0 
dx k 


( 19 ) 
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d ~ ~ — dp d ~ 

— (p W* «,) = - — + — ( % - P“k «< ) 
dx k oXi dx k 


a. io) 


d 

dx k 


(pu k h ) = 


d 

dx k 



(Ml) 


Tt mav be noted that the equation of state 


P 


N 

pRoT^ 

( = 1 


21 

W, 


becomes 

N _ 1 

P = R 0 V ( PTY. + pT"Y" ) - 
. . vv, 


(M2) 


(M3) 


When the molecular weights of different species do not vary greatly, one can write 
p=pRT (114) 


where R is the mass-based gas constant applicable to the system under consideration. 

It has now become well recognized that several difficulties arise in the course of 
Favre-averaging. Such difficulties pertain in general to the various nonlinear terms 
that arise in the describing equations. Unfortunately, there is no clear method of 
resolving the difficulties on a reasonably universal basis even with an appeal to 
experimental results. Yet another problem in Favre-averaging is that the molecular 
transport terms need to be modeled specifically for the Favre- averaged form of the 
equations, and such models cannot be the same as in the original equations. 
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2. Simplified Equations with No Variations in Density 

We consider in the following a particular case of a flowfield with temperature (or 
concentration fluctuations) in which both the flow Mach number and the density- 
variations are negligibly small. 

Equations for mean and fluctuating quantities may then be written as follows. 
Mean flow 

dUi „ dUi dp d “,u k d 2 Uj 

h Uv = h V =- 

dt dx k dxi dx k dx k 

Mean temperature 

dT ,, dT b 2 T 

— + Uk — — 1 - K » 

dt dx k dx k dx k 

Intensity of velocity fluctuations 


a. i5) 


(1.16) 


d_ 

dt 




d -r~ * ~i 

+ — [p »k + - q »k 
dx k \ 2 


dUj d\ 

Ui“k — + vu, — 
dx k dx k 


a. i7) 


Variance of temperature fluctuations 



(1.18) 
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Mean heat flux 


d_ 

~dt 


— d — d — dT — dU, n fy 

(Out) + U k — ( 0u ,) + — (Ou,u k ) = - UjU k - Ou k - 6 — 

dx k ax k ox k ox k dxj 




de „ du, 

1- vQ — 

dx k dx k 


(v + K) 


(h^dd- 
dx k dx k 


(1.19) 


It may be noted that p' represents the pressure fluctuation. Further, it may be of 
interest to observe that the second term on the right hand side of Equation (1.17) 
represents the rate of destruction of temperature fluctuation by molecular diffusion. 

Fluctuating flow 


dt dx k 


1 dp' dU, d , % 

- — u,u k ) + 

p oXj ox k axj 



( 1 . 20 ) 


Fluctuating pressure governed by Poisson equation 


d 2 p , dU t du, d 2 , 

— V = 2 — - — + (UjUj 

dxj dx k dx k dx k dxj 


UjUj) 


Fluctuating temperature 


de T , de 
— + U k — 
dt dx k 


d — d 2 T 

{ufi — ufi) + K — — j- 


dx, 


dxi 


( 1 . 21 ) 


( 1 - 22 ) 
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APPENDIX n 

APPLICATION TO MIXING LAYER 

We consider the mixing layer studied originally by J. Laufer (NACA Report 1174, 
1954) and H. Liepmann. The layer is formed by a single stream flowing over a 
stationary fluid. In the problem considered here, it is assumed that the stationary 
fluid is at a higher temperature then the flowing fluid, while both fluids are 
incompressible and at constant densities. 

The plan of attack has been as follows: (a) to obtain the mean velocity field and 
the shear stress distribution in the mixing layer; (b) to invoke spectral analogy 
between velocity and temperature fields and thereby establish the mean temperature 
field and temperature variance distribution; and finally, to predict the velocity- 
temperature correlation distribution which could be checked against experimental 
results. It may be said at the outlet that no calculations could be completed within 
the project period. 

H.l. Formulation of the Problem. 

Referring to Appendix I, Eqs. 1.15, 1.20 and 1.21 can be solved to obtain mean 
flow, turbulence intensity of velocity of fluctuations and Reynolds stress. The 
modelling required pertains to eddy-eddy interactions. Eq. 1.17 can be utilized to 
check the modelling. 

In Eq. 1.18, the analogy between velocity and temperature intensities can be 
invoked. One then solves Eqs. 1.18, 1.22 and 1.16 to obtain ur(? and mean 
temperature as functions of xj. Returning to Eq. 1.22 it is possible to establish the 
skewness factor and turbulent conductivity parameter that are compatible with the 
analogy between velocity and temperature. 

Finally, it is possible to extract the pdf of velocity and temperature fluctuations 
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and compare them with measured data. 

In analyzing the flowfield one may start with assuming that the flowfield is 
homogeneous in the spanwise direction and then proceed to a formulation in which 
the flowfield is fully inhomogeneous. Such a procedure is especially useful in 
extracting physical structures. 

It may be pointed out that the entire formulation is based upon concentrating 
attention on the largest eddy represented by the first mode and a few wave numbers 
in the overall spectrum. Therefore, the analogy between temperature and velocity 
intensities may in fact be partially valid. This may cause errors in the individual 
contributions of production, dissipation and transport to the stresses and the 
velocity-temperature correlations. A multi-modal analysis, and on the other hand, 
requires some experimental input. 

n.2. Orthogonal Decomposition and Spectral Analysis 

The orthogonal decomposition and spectral analysis can be illustrated utilizing the 
case of velocity. As stated elsewhere, the instantaneous velocity can be written in 
terms of a time-mean value and its fluctuation: 

Uj = Uj 4 - Uj (II. 1) 

where the velocity fluctuation can be represented in a series in terms of orthogonal 
functions: 


U,(x.[) = £ a n (t) <*>/'" (j) n = 1. 2. 3 (II.2) 

n = i 


The random coefficients, a n (t), are uncorrelated with each other. Hence, 
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a m a„ = X 2 «) • 6 mn , b mn = kronecker delta, 


OI. 3 ) 


the overbar representing an ensemble average. Assuming that the determinate 
are orthonormal functions 


J *<'W'dx = 5 „ 01 . 4 ) 

It then follows that the two-point velocity correlation, Rjj, can be expressed in the 
form, namely: 


R,j (x,x';t) = £ X 2 (x) • 4>l n) (*' ) 01.5) 

n ~ 1 

Now, one can set up dynamical equations for the n tl1 mode by an appeal to 
Navier-Stokes equations. After some manipulation, one obtains the following 
dynamical equation. 


djy , „ ,9U, ( " " 

dt J dx, dxj { ", “ 

J J ^ p-\ q= I 








dp a 2 $7" ) 

— h V =— 

Ax, dxj 


01.6) 


where 

$l n) = VX^> and = 

p Vx<"> 

Also, the Poisson equation for pressure fluctuation becomes the the following. 
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3 2 7T ^ _ 2 3 Uj 3<j> % * f y^ y^ Qf i a /> a <y ^(fl) 0^) 

dxf dx k dxj dx k d Xj [x^X^X^]* 7 * . 


The continuity conditions yields 

°±L - -Q (II. 8) 

dxj 

Eqs. n.e-n.8 can be specialized to the first mode, n = 1, considered on the large 
eddy. 

The eddy-eddy interactions are then modelled as follows. 


E E 


p = i ?=> 


ot | 0 <p) 

X, 



A (l) . (l) 
<Pi <f>j 



(II -9) 


where 



X, 

the skewness factor of the random coefficient; and is a form of damping coefficient 
referred to as eddy viscosity for convenience. 

II. 1.1. An Alternative Formulation 

At this stage, we have also considered an alternative formulation. It consists in 
including two types of eddies, one that may be referred to as large and the other as 
small, it is assumed that the large and the small eddies are sufficiently separated. 
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The velocity fluctuation is then written in the following form. 


u,(x,t) = a<t>i (x,t) + a'<P,'(x,t) (11.11) 

The two deterministic functions are orthogonal in the (x,t) space of each other. The 
random coefficients do not have any correlations whatsoever between them, that is 


aa’ = 0 

The Reynolds stress tensor then becomes 


u,Uj — a 2 <t>j 4>j + a' 2 4>i 4>j— (Vo 5 </>,) (a' 2 0 ; ) 


( 1112 ) 


The approximation is an important one to proceed further. 

The eddy-eddy interactions again involve a skewness factor and an eddy viscosity. 
The second-order eddy viscosity tensor ey is written in the following form. Defining 


5 , = 


a 2 a' 


a 2 y U ' 2 


So = 


aa 


a’ 2 'To 2 


we write 


$,(<*>' <j>j ~ <t>i <t>j) - S 2 (4>; <!>’) = - e ik 



01 . 13 ) 


There are some complex arguments in setting up Eq. (11.13). In any case, this 
approach has not been developed further. 
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n.1.2. Temperature Fluctuations 

In a manner similar to the foregoing, it is also possible to decompose the 
temperature fluctuation, 9 and the correlations Uj$. Thus 

oo 

0 (*/)=£ 0„(f) 7 <n> « = 1,2,3,.... (11.14) 

n = 1 


and 


«,* = £ o n ^ n) (x)y^ (11.15) 

n — ! 

A dynamical equation can also be constructed for 9 corresponding to Eq. II. 6. 
However, as stated earlier, it is proposed to invoke a similarity between temperature 
and velocity intensities. 

n.1.2. Spectral Decomposition 

Assuming that the flowfield is homogeneous in the z-direction, spanwise with 
respect to the developing mixing layer, one can construct the spectrum as follows. It 
may be observed that a more general formulation in which the flowfield is 
inhomogeneous in all three directions is easily set up formally. However, it is 
worthwhile starting with the simpler formulation. 


4>, (I) (* 3 .‘ x.y.z.t) 


1 

2ir 


0, (l) (x.y.z.t) exp (-ik 3 z) dz 


(11.16) 


where k 3 is the (scalar) wave number corresponding to the z-direction. For any z, 
then, turbulence contains various length scales whose sizes are proportional to l/k e . 
The spectrum, <f>{^ (k 3 ), is the integrated value of the first mode, (z), 
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corresponding only to l/k 3 , over the entire homogeneous field, z. Inversely, one can 
write 


& (,) (x,y,z.t) 



4>- l) {k 3 ; x.y.z.t) exp (ifc 3 z) dz 


01 . 17 ) 


The spectrum is a complex variable which may be divided into real and imaginary 
parts as follows. 


i, = (<£,)* + Hi,) i 


Proceeding similarly with respect to fr and 5, one can write 


~ P[ + iP 2 

0 2 = 7*3 + ‘P 4 

03 ~ 7*5 + lP 6 
■K = P-i + fPg 

e=p 9 + fp w Oils) 

It is important to recognize here that, while the flowfield is assumed to be two- 
dimensional, with homogeneity being imposed in the third direction, turbulence is 
treated correctly as being three-dimensional. 

Various quantities such as UjUj, uf, q 2 , uv/q 2 , orientation of principal axes and 
anisotropy can be expressed in the mixed (x,y, k 3 ; t) space. Similarly, the various 

turbulence quantities related to temperature, 8 2 , u,0, can also be expressed in spectral 


space. 
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n.2. Analogy Between Temperature and Velocity 
The analogy is expressed in the following form 




( 11 . 19 ) 


where j3 is in general of 0(1), but is taken as equal to unity in the first instance. 

As stated earlier, q 2 is expressed in spectral form in the first instance, while 
simultaneously establishing the mean flow distribution. The spectrum of e 2 can be 
related to the spectrum of q 2 but only with a knowledge of cFI/dy. thus, it is 
necessary to solve Eqs. 1.18, 1.22 and 1.16 simultaneously to determine mean 
temperature and u^#. It may be pointed out that in accomplishing this in practice, it 
is useful to equate <9T/dy to dU/dy and thus use the spectrum of q 2 as the spectrum 
of 9 2 . One can then proceed to modify the analogy through successive iterations. 
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